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Abstract: An extreme ultraviolet solar corona multispectral imager can allow direct observation of high temperature 
coronal plasma, which is related to solar flares, coronal mass ejections and other significant coronal activities. This 
manuscript proposes a novel end-to-end computational design method for an extreme ultraviolet (EUV) solar corona 
multispectral imager operating at wavelengths near 100 nm, including a stray light suppression design and computational 
image recovery. To suppress the strong stray light from the solar disk, an outer opto-mechanical structure is designed to 
protect the imaging component of the system. Considering the low reflectivity (less than 70%) and strong-scattering 
(roughness) of existing extreme ultraviolet optical elements, the imaging component comprises only a primary mirror and a 
curved grating. A Lyot aperture is used to further suppress any residual stray light. Finally, a deep learning computational 
imaging method is used to correct the individual multi-wavelength images from the original recorded multi-slit data. In 
results and data, this can achieve a far-field angular resolution below 7", and spectral resolution below 0.05 nm. The field of 


view is +3 Ro along the multi-slit moving direction, where Ro represents the radius of the solar disk. The ratio of the 
corona's stray light intensity to the solar center's irradiation intensity is less than 10~ at the circle of 1.3 Ro. 
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1. INTRODUCTION 


Solar coronal observation is very important for the protec- 
tion of both spacecraft and planet Earth. The EUV band is sig- 
nificantly different to visible light and other bands of electro- 
magnetic radiation, with a shorter wavelength and higher pho- 
ton energy!!-3]. Because this band is strongly absorbed by 
Earth’s atmosphere, EUV telescope imagers work best in 
space. An EUV solar corona imager can enable direct observa- 
tions of high temperature coronal plasma, and the most attrac- 
tive for this application are multispectral imagers, consisting 
of a wide field imager and a spectrometer. Most materials 
absorb EUV radiation, preventing the use of prisms and trans- 
missive optical elements. Additionally, due to the low EUV 
reflectivity of materials, it is unrealistic to add too many 
reflective elements!*+ >] into EUV multispectral imagers. For 
example, a basic reflective Fourier transform imaging spec- 


trometer, in which light is reflected several times, is inappro- 
priate for EUV lightl*8], Some EUV thin-film filters have 
been used in the EUV imager on the Solar Dynamics Observa- 
tory satellite and the Atmospheric Imaging Assembly satel- 
lite, but these only have a spectral resolution of approximately 
0.5 nml?!21, 

With advancements in multilayer film and grating 
fabrication techniques, curved gratings have been introduced 
into EUV spectrometers. These curved gratings have unique 
diffraction capabilities, enabling spectral imaging and disper- 
sion with only a single grating, reducing the system's reflec- 
tive and diffractive surfaces significantlyl!3-15], Conse- 
quently, an EUV solar corona multispectral imager would 
combine a telescope system with a curved grating!!*!8], The 
telescope images the incident light into a slit, which serves 
as both the field stop for the telescope and the slit for the spec- 
trometer!!°], The curved grating here works simultaneously 


as both imager and grating, removing the requirement to colli- 
mate mirrors, thus greatly improving the system sensitivity. 
The system can achieve a wide field-of-view (FOV) 
through scanning movement of the slits!?-27], The 
CORONA project (America’s first satellite program) 
employed a planar grating combined with a focusing 
mirror??3-25], able to perform multispectral imaging with 
high spectral resolution and wide FOVI?*. 

For an EUV solar corona multispectral imager with 
high sensitivity and efficiency, we propose an end-to-end 
computational design method with separate optical elements 
and a deep learning computational image recovery algo- 
rithm. For suppressing stray light from the solar disk, the 
external occulter (EO) is innovatively designed, with an annu- 
lar light aperture and an inclined spherical mirror with a cen- 
tral opening. For the imaging component, the simplest config- 
uration is adopted, with only a primary mirror and a curved 
grating, which can achieve both high spatial and spectral reso- 
lution in a wide FOV. A unique multi-slit arrangement 
greatly improves the scanning efficiency. In addition, 
enhancement by deep learning (DL) methods can correct dis- 
tortion caused by the curved grating. The entire system com- 
prises only two optical reflective surfaces, ensuring high trans- 
mission efficiency. 


2. RESULTS 


2.1. The EUV Solar Corona Multispectral Imager 
System 


The main specifications of the EUV solar corona multi- 
spectral imager system are listed in Table 1. The EUV spectral 
ranges is 94 nm to 109 nm, namely 94 nm, 97.5 nm, 
97.7 nm, 102.5 nm, 103.2 nm, 103.7 nm and 109 nm. The 
half-angle of FOV along the multi-slit moving direction is 
3Ro and the half-FOV along the slit is 2 Ro. The spectral res- 
olution is better than 0.05 nm, and the far-field angle resolu- 
tion is better than 7". The design uses 5 slits for parallel-scan- 
ning spectral imaging. The design of the EUV solar corona 
multispectral imager is shown in Fig. 1, including an EO 
structure, entrance pupil stop, off-axis primary mirror, multi- 
slit plane, spectral dispersion component, and Complemen- 
tary Metal-Oxide-Semiconductor (CMOS) image sensor. 
As stray light beams from the solar disk are very strong, 
they should be prevented from entering the aperture, so an 


outer stray light suppression component is required before 
the imaging system. The primary mirror is then placed after 
the entrance pupil stop. The primary mirror is an off-axis 
parabolic mirror and is used to capture the first image of the 
solar corona, with the 5-slit multi-slit mask placed at the imag- 
ing plane of the primary mirror. A curved grating is used 
for multispectral imaging, conjugately imaging the first 
image plane to the CMOS image sensor. A Lyot stop is 
placed at the exit pupil to suppress stray light. The CMOS 
image sensor is placed on a tilted plane. Fig. 2 shows the 
quantitative imaging performance of the system, including 
root mean square (RMS) data, spot radii, and modulation 
transfer functions (MTFs). In Fig. 2A, B and C, the RMS 
spot radii are all below 10 um, where the pixel size of the 
CMOS image sensor is 25 um. Additionally, the MTF is better 
than 0.2 at 40 lp/mm across all FOVs, at wavelengths of 94 
nm, 102.5 nm, and 109 nm. 


2.2. Stray Light Suppression 


To prevent stray light radiation from the solar disk, an 
outer stray light suppression shielding system is employed. 
The opto-mechanical structure is shown in Fig. 3A, with 
two separated components. The front body includes a circular 
obscuration aperture (COA), while the rear body employs 
an annular spherical mirror aperture (ASMA). Fig. 3B 
shows the irradiance distribution of the annular incident win- 
dow, and Fig. 3C shows the irradiance distribution on the pri- 
mary mirror of the coronagraph, within the FOV of 0-2 Ro. 
Results indicate that rays within the 0-1.3 Ro FOV are 
imaged at the annular aperture and kept out of the imaging 
apparatus, where the ratio of the corona's stray light intensity 
to the solar center's irradiation intensity is less than 10~6 at 
the circle of 1.3 Ro. 

There are two main sources through which stray light 
can be introduced into the imaging system, which can 
reduce the image contrast of the CMOS image sensor. One 
source is light diffracting and scattering at the edge of each 
element. Another, as the EUV wavelength is short, is from 
light scattering off any roughness or defects in the mirror 
and grating. To minimize the effects of this, a Lyot stop is 
placed at the exit pupil of the imaging component, as shown 
in Fig. 4A. Calculations show that the diffracted light forms 
an annular illumination pattern with a diameter of 8.9 mm 
at the location of the Lyot stop. Consequently, the optical 
diameter of the Lyot stop is set to 8.7 mm to completely 


Table 1. EUV solar corona multispectral imager system specifications 


Wavelengths 


Entrance pupil 
Obscuration half-FOV 
FOV in multi-slit scanning direction 
FOV along the slit 
Focal length 

Spectral resolution 
Spatial resolution 

Number of slits 
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94 nm & 97.5 nm & 97.7 nm & 102.5 nm & 
103.2 nm & 103.7 nm & 109 nm 


50 mm 
1.3 Ro 
+3 Ro 
+2 Ro 
755 mm 
<0.05 nm 
<7" 
5 
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Fig. 1. Schematic of the EUV solar corona multispectral imager. (A) Opto-mechanical structure, from the outer stray light 
suppression component to the CMOS image sensor. (B) Slit mask for multi-slit scanning used in EUV multispectral imaging. (C) Optical 
light path, including outer stray light suppression and imaging components. 
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Fig. 2. Image performance of the EUV solar corona multispectral imager system. (A) RMS spot radii change with wavelength under 
different FOVs. (B)-(C) RMS spot radii versus FOV at different wavelengths. (D)-(F) MTF of the optical system at wavelengths of 94 


nm, 102.5 nm, and 109 nm. 


block the diffracted light. The Lyot stop is positioned 
178.2 mm before the image plane with a tilt angle of —6.2°. 
The system vignetting caused by the Lyot stop is shown in 
Fig. 4B. The maximum unvignetted light fraction of the sys- 
tem is over 95%, which will not influence the image qual- 
ity. Finally, stray light is suppressed less than 10~ at the circle 
of 1.3 Ro. 


2.3. Primary Mirror 


A primary mirror is used for the first imaging, which is 
also the multi-slit plane. Here, an off-axis parabolic mirror 
works as the primary mirror, similarly to a telescope objec- 
tive, the design parameters of which are given in Table 2. 
The radius of the parabolic mirror is —8 10 mm, the conic coef- 
ficient is —1, the decentration is 50 mm, and the clear aperture 
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Fig. 3. The EO. (A) Mechanical structure of optics. (B) Irradiance distribution at the annular aperture. (C) Irradiance distribution at the 
image plane of the primary mirror and the multi-slit plane. 
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Fig. 4. Lyot stop to suppress stray light rays in the imaging 
system. A: Position of the Lyot stop. B: Vignetting calculation 


of the system caused by the Lyot stop. Please note that the y- 


axis here shows only 0.9-1. 
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Fig. 5. Sagittal height of the primary mirror. (A) Surface sagittal height diagram. (B) Cross-sectional sagittal height diagram. 
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Position/mm 


Table 2. Primary mirror parameters 


Radius of Conic : Clear 

. Decentration 
curvature coefficient aperture 
—810 mm =f 50 mm 55.4 mm 


is 55.4 mm. The surface sagittal height profile and the cross- 
sectional sagittal height profile of the primary mirror are 


shown in Fig. 5. 
2.4. Curved Grating 


The parameters of the curved grating are divided into 
two parts: surface shape parameters and engraved line parame- 
ters. The surface shape is composed of an ellipsoid base sur- 
face superposed with XY extended polynomials. The surface 


shape expression is expressed as 


cu? a 
"A e e, +) AE; (x,y), 
1+ Vl- u2 iÍ 
where, u can be expressed as 
w =a x+ b’y’. 
Boo 
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Position/mm 
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(2) 


In Eq. (1) and Eq. (2), x, y and z are the intercept coordi- 
nates of the intersection point between the light ray and the 
surface, while a, b, and c are the coefficients defining the 
major and minor axes of the elliptical shape. The grating 
line density is given by 

1 


1 
= =de = — +y +py +yy’. 


T To ©) 


The specific surface parameters of the grating are 


shown in Table 3, Table 4, and Table 5. The surface sagittal 
height diagram and the cross-sectional sagittal height diagram 


of the grating are shown in Figure 6. 
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Table 3. Grating surface parameters 
x tilt 
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Table 4. Surface sagittal height polynomial parameters 


X1Y0 5.208x 10-6 X1Y1 —3.318x107 
X0Y1 —8.697x1076 X0Y2 11.416 
X2Y0 3.2 X3Y0 1.245x10~ 
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Fig. 6. Sagittal height diagram of curved grating. (A) Surface sagittal height diagram. (B) Cross-sectional sagittal height diagram. 


Fig. 7 displays the grating line density corresponding to 
the center position of the spot, indicating that the groove den- 
sity varies with the relative position in the Y direction. The 
relative landing point position of the edge ray of the target 
FOV on the diffraction grating surface is —13.31 mm. The cor- 
responding grating line density at this point is 2034 
lines/mm. The relative landing point position of the lower 
edge ray of the target FOV on the diffraction grating surface 
is 13.31 mm, the grating line density corresponding to this 
point is 1970 lines/mm. Using Eq. (3) to integrate the spot 
range, the total number of slots obtained is 36 432, which 
meets the system's requirements for spectral resolution. 


2.5. Computational Multispectral Image Correction 


Each original spectral image is stitched from slit-bar 
data, which is directly recorded by the CMOS image sensor. 
These stitched images contain distortion commonly found 
in slit-scanning spectral imagers. A DL method with a convo- 
lutional neural network (CNN) is used to correct imaging dis- 
tortion for the stitched spectral images, as shown in Fig. 8. 
We use a total of 5000 sets of data to train and test the perfor- 
mance of our computational multispectral image recovery, 
constructed within the PyTorch DL framework. The input 
and target values of the CNN are distorted corona images 
and ideal coronal images respectively. 4500 sets are used 
for training the CNN network and 500 sets are used for test- 
ing. A total of 20 epochs are performed, and the entire calcula- 
tion process is performed on a desktop computer running 
the Microsoft Windows 10 operating system, with Core i7- 
7700K CPU @ 4.2GHz (Intel) with 64GB RAM and 
GeForce RTX2080Ti GPU (NVIDIA) hardware. The training 


time is approximately 2 hours, and the correction time of a 
single image is less than 200 ms. Test results are shown in 
Fig. 8. In 500 sets of test data, the average structural similarity 
value of the corrected images and raw image reaches 0.898, 
which is higher than the 0.769 of the distorted images. 
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Fig. 7. Grating line density distribution map. The diameter of 

the curved grating is 40 mm. As examples, the point position at 

—13.31 mm on the diffraction grating surface has a grating line 

density of 2034 lines/mm, while the point position at 13.31 mm 

has a grating line density of 1970 lines/mm. 


3. DISCUSSION 


There are two important considerations for an EUV 
solar corona multispectral imager, which can easily be over- 
looked. The first is to suppress stray light. Light from the 
solar disk is very strong, and this must be excluded from the 
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Y axis 


——Ea 


Fig. 8. Deep learning distortion correction for multispectral images correction. The reddish-brown area in the figure shows 
perfectly circular masks. (A) Imaging FOV diagram. (B) Simulated solar corona raw image. (C) Stitched image after slit scanning, 
from slit-bar data at 94 nm, with distortion. (D) Output of the deep learning CNN, without distortion. 


imaging component. Additionally, to avoid internal scattering 
within the apparatus, the edges of all apertures should be 
crafted with minimal imperfections. All mechanical struc- 
tures should be coated in light-absorbing material to minimize 
reflected light. Meanwhile, the mechanical structure should 
be strong and any vibrations should be sufficiently dampened 
to survive a rocket launch. Therefore, the mechanical struc- 
ture should be designed with multiple chambers, with good 
inner space distribution, to give the instrument increased sta- 
bility. This can also limit the propagation of scattered light 
within each narrow individual chamber, preventing the prolif- 
eration of stray light. 

The second consideration is to design and arrange suit- 
able slits, as detailed in Section 5. The slit width should 
match the optical resolution and the pixel size of the CMOS 
image sensor. If the slits are too narrow, diffraction and expo- 
sure time can become problematic, but slits which are too 
wide will reduce spatial resolution to unacceptable levels. 
Additionally, the energy distribution of the EUV spectrum 
peaks at a specific wavelength. The arrangement of slits 
must be considered according the wavelength distribution pro- 
duced by the curved grating. 


4. CONCLUSION 


In this paper, we propose an end-to-end computational 
design method for an EUV solar corona multispectral 
imager with high sensitivity and efficiency. It uses a design 
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with separate optical elements and a DL computational 
image correction algorithm. The design uses an external EO 
with an annular light aperture and an inclined spherical mirror 
with a central opening, to suppress stray light from the solar 
disk, with a unique multi-slit arrangement to greatly 
improve scanning efficiency. The entire system comprises 
only two optical reflective surfaces, ensuring high transmis- 
sion efficiency. Following image acquisition, distortion 
from the curved grating can be corrected using a CNN. 
From the design results, the EO enhances stray light suppres- 
sion, and the spectrometer exhibits outstanding spatial and 
spectral imaging performance. The far-field angular resolu- 
tion of this design is better than 7", and the spectral resolution 
is better than 0.05 nm. The FOV is +3 Ro along the multi- 
slit moving direction. The ratio of the corona's stray light 
intensity to the solar center's irradiation intensity is less than 
10% at the circle of 1.3 Ro. We believe that this novel 
method can lead to the design and development of further 
new EUV spectral imaging systems. 


5. METHODS 


5.1. End-to-end Computational Design 


An end-to-end computational design for the EUV solar 
corona multispectral imager is presented in Fig. 9. The physi- 
cal imaging model is generalized, from the object ‘solar 
corona’ to the final recovered images. In Fig. 9B, a stepwise 


Deep learning CNN 


» 


B 


External occulter 


Telescope imaging system 


Multi-slit offieimaseiplane 


Spectral dispersion. 
system 


Fig. 9. End-to-end computational design for the EUV solar corona multispectral imager. (A) The end-to-end roadmap from the 
object ‘solar corona’ to the final recovered images. (B) Each design point, including stray light suppression, telescope objective, multi- 
slit mask, curved grating and the deep learning recovery algorithm, along the end-to-end roadmap. 


approach is adopted based on imaging requirements, culminat- 
ing in the final imaging criteria. To block the solar surface 
radiation from the coronagraph system as early as possible, 
the design of the EO and the vignette analysis are carried 
out first. Then, to obtain a clear image of the solar corona at 
the slit position, the telescope primary mirror is designed com- 
putationally. To reduce the full-field scanning imaging time 
and eliminate the overlap of imaging spectral lines in different 
FOV, it is necessary to determine the number of slits and 
slit arrangement. Finally, DL computational image recovery 
is used to correct aberrations and improve image quality. 


5.2. Design of Stray Light Suppression Components 


A schematic diagram of the EO stray light suppression 
structure is shown in Fig. 10. The structure consists of a 
COA and an ASMA with an inclination angle y. The center 
of the COA is an occluder, and the center of the ASMA is 
an aperture with a radius of rọ. A spherical mirror with a spe- 
cific tilt angle, y, can image the sun at the annular optical 
entrance, thereby reflecting it out of the solar corona 
imager. The inner radius of the COA is rı, the outer radius 
is r2, and the positions Eg and Mọ are separated by dis- 
tance, L. If the image of the sun is located on Eo, the distance 


Image of 
the sun 


Annular spherical 
mirror aperture (Mo) 


Circular obscuration 
aperture (Ep) 


Fig. 10. EO stray light suppression structure and design. It 
consists of a COA and an ASMA with an inclination angle. 
The center of the COA is an occluder, with an aperture at the 
center of the ASMA. The spherical mirror images the sun at 
the annular optical entrance, thereby reflecting it out of the 
solar corona imager. 


d between the image center and the COA center needs to sat- 
isfy the formula 


rth<d< rm-h, (4) 


where, A is the half-height of the image of the sun, so that ¥ 
satisfies the formula 
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1 1 
5 aretan[(r1 +h)/L]< Y< qarctan| (72 —h)/L). (5) 
Here, 9 is the half-FOV angle of the solar surface 
image, and @; is the maximum half-FOV angle of the solar 
coronal image. All light rays from 0-69 will be intercepted 
by the ASMA, will be imaged on the surface of the annular 
entrance where the COA is located, and will be reflected 
out of the solar corona imager. All the solar coronal light 
from 9-6, enters the system, while the light greater than 0; 
will be partially intercepted by the ASMA. If the pointing 
accuracy of the satellite is g, then 7; and r2 satisfy the equa- 
tion 
rı = ro + LX (tan 6o + tan y) 7 
r > ro + Lx (tan, +tang) ` (6) 


5.3. Primary Mirror Design 


To meet the desired spatial and spectral resolution 
requirements, clear imaging at the primary image plane of 
the system is needed, i.e., at the slit position. The main mirror 
adopts an off-axis parabolic design!?”]. When determining 
the parameters of the off-axis paraboloid, the following fac- 
tors should be comprehensively considered: (1) Since a 
multi-slit plane is placed at the primary image plane of the pri- 
mary mirror and an image of the corona region is obtained 
through scanning the multi-slit device, the imaging focusing 
performance (spot focusing diameter) of the primary mirror 
needs to match the width of the slit. (2) The focal length of 
the spectral imaging system only depends on the focal 
length of the primary mirror and the magnification of the 
curved grating. Therefore, the focal length design of the off- 
axis parabolic primary mirror needs to match the spatial reso- 
lution!28], (3) To avoid mechanical interference between the 


optical path and optical elements, there should be sufficient 
space for mechanical structure assembly and adjustment. 
The off-axis parabolic mirror needs to be sufficiently off- 
axis. At the same time, the off-axis paraboloid serves as the 
front telescope system of the spectral imager, and its design 
needs to meet the imaging focusing performance of the target 
FOV at the final image planel??-3!1, 


5.4. Slit Distribution Design 


The multi-slit plane is located at the primary image 
plane and consists of straight slits with unequal spacing, as 
shown in Fig. 11. From the focal length of the telescope imag- 
ing system, the field diaphragm size at the first image plane 
can be calculated as 

{ Yrs = 2 fr Xx tan wy (7) 
Xrs = 2 fr X tan wx 


where, Yps is the width in the Y direction and Xps is the 
length in the X direction. In Fig. 11, Nw+(N-1)ds; = 
Yrs, where w is the slit width, d,; is the slit spacing, and N 
is the number of slits. The image of the slit through the grating 
has a width wı on the CMOS image sensor. The spectra pro- 
duced by w; will partially overlap, limiting the resolvable 
spectral bandwidth. The spectral resolution corresponding 
to the geometric width of the slit image should satisfy 


do cos @cos y 
| ES 


da=w (8) 


mrg 

Multiple slits enable different spectral line segments at 
different FOVs to be imaged on the CMOS image sensor 
simultaneously. To reduce the time cost, multiple slits 
should be involved in imaging while ensuring spectral resolu- 
tion. The slit spacing ds ; can be calculated with the formula 
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Fig. 11. Schematic diagram of the multi-slit mask, where ‘px’ denotes size in pixels. The spectral dispersion distributions 
corresponding to the five slits are independent and are located together on the CMOS image sensor. The spectral bandwidths 
will partially overlap, but the spectral line positions do not overlap. 
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dsi =Bdmji i=1,2,---N, (9) 


where, ds, is the spacing between i slits, dm,;is the minimum 
spacing between slits i and i+1 on the image plane, and £ is 
the magnification of the curved grating. 


5.5. Curved Grating Design 


To correct chromatic aberration and obtain high spectral 
resolution imaging, the spectral dispersion system uses a vari- 
able line spacing curved grating. The aberration correc 
tion principle diagram is shown in Fig. 12. The point 
Ao (Xa, Ya,0) is the center of the slit, and AAo is any off-axis 
light source parallel to the grating line. The point P (x,y,z) 
is any point on the grating surface. After the incident light 
rays AP and AO pass through the curved grating, they are 
imaged at point B (xp, Yb,Zb). Bo(Xp,¥p,0) is the projection 
of point B on the plane XOY. The angle i between the light 
AoO and the X-axis is the incident angle, and |OAo| = ra is 
the incident arm length of the grating. The angle 0 between 
the light BoO and the X-axis is the diffraction angle, and 
|OBo| = rg is the exit arm length of the grating. The surface 
sag of the grating can be expressed using Eq. (1). The density 
of the grating changes along the Y-axis. The total number of 
gratings can be expressed as 


¢ ba 4 


2,6 3 


(10) 


NO) =(y+ Sy + pay tR ++} do, 


=> 


2 
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Fig. 12. Diagram of the curved surface variable line pitch 
grating. 

where, N (y) is the total number of score lines accumulated 
from the origin O to any point P in the Y direction, dois the 
score line spacing at the origin O, and @; is the spatial variation 
parameter of the score line density. Bringing the coordinates 
of each point into the optical path function of APB and 
expanding it into the Taylor series of y and z, we obtain 


F(x,y) =ra +rg+ Y yietFy 
Tk . (11) 
Fik = C jx + (mA/do) M jx 


In Eq. (11), j and k cannot simultaneously be zero, and 
each term corresponds to a grating aberration of the curved 
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Fig. 13. Schematic diagram of our deep learning CNN. The input is a distorted stitched solar coronal image and the output is a 
corrected solar coronal image. The network layers used in the CNN structure are shown to the right. 
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grating. By optimizing the spatial variation parameters of 
the curved grating, the spectral focus curve and the spatial 
focus curve can be made to intersect at the astigmatic point, 
thereby achieving non-atigmatism imaging. Through the opti- 
mization of superposition polynomials, high-order aberra- 
tions can be corrected and the optimal solution of the system 
can be obtained. 


5.6. Computational Imaging Method 


The stitched image will inevitably have distortion, 
caused by the slit-bar data recorded with the CMOS image 
sensor, like a smile curve. Although small, this distortion is 
still non-negligible. This is caused by the inherent characteris- 
tics of the curved grating and off-axis imaging. Therefore, it 
is necessary to perform computational post-processing on 
the distorted stitched solar coronal image to obtain a corrected 
solar coronal image. The essence of the image distortion cor- 
rection problem is to spatially transform the distorted 
image, and the transformation relationship is fixed in this sys- 
tem. A DL CNN can be used to fit this transformation!32-34], 
which is shown in Fig. 13. The entire process can be viewed 
as feature extraction of the input image and pixel-by-pixel 
regression of the output results!?>- 3°], Therefore, the CNN 
first convolves and down-samples the input data, performs 
feature transformation in 9 residual blocks, and finally up- 
samples through transposed convolution to reconstruct a dis- 
tortion-free image. When training the model, we use the 
mean square error loss as the loss function, which is 
expressed as 


n 


1 Fey) 
MSE=- X (Yi-?), 


i=1 


(12) 


where, Y; is the ideal coronal image, Y; is the reconstructed 
coronal image, and n is the number of elements of the output 
tensor. 


6. AVAILABILITY OF 
MATERIALS 
All data needed to evaluate the conclusions in this 


study are presented in this manuscript. Additional data 
related to this paper may be requested from the authors. 


DATA AND 


7. ABBREVIATIONS 


ASMA: Annular Spherical Mirror Aperture 

CME: Coronal Mass Ejection 

CMOS: Complementary Metal-Oxide-Semiconductor 
CNN: Convolutional Neural Network 

COA: Circular Obscuration Aperture 

DL: Deep Learning 

EO: External Occulter 

EUV: Extreme Ultraviolet 

FOV: Field Of View 

MTF: Modulation Transfer Function 
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RMS: Root Mean Square 
SDO: Solar Dynamics Observatory 
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